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Qh| We discuss the effect of the instanton induced, six-fermion effective Lagrangian on the decays of the lightest 



scalar mesons in the diquark-antidiquark picture. This addition allows for a remarkably good description of light 
scalar meson decays. The same effective Lagrangian produces a mixing of the lightest scalars with the positive 
parity qq states. Comparing with previous work where the qq mesons arc identified with the nonet at 1200-1700 
MeV, we find that the mixing required to fit the mass spectrum is in good agreement with the instanton coupling 
obtained from light scalar decays. A coherent picture of scalar mesons as a mixture of tetraquark states (dominating 
qq \ in the lightest mesons) and heavy qq states (dominating in the heavier mesons) emerges. 
(N| ' PACS 12.38.Aw, 12.39.Mk, 14.40.-n 
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3 ■ 1 Introduction 

00 ' 

We study in this paper the strong decays of the lightest scalar mesons: a, k, fo, ao and the relations of the light 
scalars to the scalar mesons observed in the 1-2 GeV range. 

Recent experimental and theoretical evidence for the existence of a and k [1-4] (see also [5]) indicates that light 
scalars make a full SU(3) flavor nonet. Their mass spectrum, with the peculiar inversion of the k and fo or ao mass 
ordering, speaks however against the naive qq picture. The most natural explanation for such complete multiplct 
with inverted mass spectrum is that these mesons are diquark-antidiquark bound states. The KK molecular 
constitution [6], advocated to explain the degeneracy of /o/ao with the KK threshold, would lead most likely to 
incomplete multiplets. 

The picture where the light scalar mesons are diquark-antidiquark states bound by color forces has been 
discussed by several authors [7-9]. In this picture the diquarks, which we will indicate with [01O2], are in color 
3, spin S = and flavor 3, and antidiquarks in the conjugate representations. Diquark-antidiquark bound states 
(tetraquarks, for short) naturally reproduce the SU(3) nonet structure with the correct mass ordering, as indicated 
by the explicit quark composition: 

<j [a] = [ud][ud) 

k = [su][ud]; [sd][ud] (+ conjugate doublet) 
[0] _ [su] [su] + [sd] [sd] 
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While the mass spectrum of the light scalar mesons is well understood in terms of diquark-antidiquark bound 
states, the overall picture of scalar mesons is still affected by two drawbacks. First, the strong decays into two pseu- 
doscalar mesons (S — > PP) have so far escaped a satisfactory theoretical understanding in the quark rearrangement 
picture. In particular, the fo — -> 7T7t coupling is too small compared to experiments (according to the ideal-mixing 
decomposition in |T]) it should vanish) and the ao — > r/7r coupling largely exceeds its experimental value. Second, 
the identification of the qq scalar sates is an open issue. The latter must exist, as indicated by the well identified 
nonets corresponding tqj orbital angular momentum L = 1 and quantum numbers J PC = 1 ++ , l" 1 , 2 ++ [10]: all 
predicted states are unambiguously identified but for J PC = ++ . 

In this paper we show that these two main problems are solved by the instanton induced effective six-fermion 
Lagrangian [11], the same effective interaction which solves the problem of the rj—rj' masses [12,13]. Such Lagrangian 
has two important effects in scalar mesons dynamics: (i) it generates a mixing between tctraquarks and qq states, 
(ii) it provides an additional amplitude which brings the strong decays of the light scalars in good agreement with 
data. Former studies of instanton-induced effects in scalar meson dynamics can be found in the literature [14] , but 
these two effects have not been discussed before. 

The two effects induced by the effective instanton Lagrangian are closely connected. The tetraquark-qq mixing 
makes it possible to identify the heavy scalars around 1.5 GeV as predominantly qq states with a non-negligible 
tctraquark component. The latter is essential to explain the anomalous mass spectrum of such mesons, as originally 
proposed in [8]. On the other hand, integrating out the heavy qq components, the tctraquark-qq mixing manifests 
itself into the non-standard S — » PP decay amplitude for the light scalar mesons which improves substantially 
the agreement with data. These two independent phenomena lead to obtain two independent phenomenological 
determinations of the non-perturbative parameter which controls the matrix elements of the instanton Lagrangian 
in the scalar sector. The two determinations turn out to agree, reinforcing the overall consistency of the picture. 

The paper is organized as follows: in Section [2] we illustrate the two main effects of the instanton interaction 
in the scalar sector, constructing the corresponding effective Lagrangians in terms of meson fields. In Section [3] 
we present a numerical analysis of S — PP decays in this scheme, demonstrating the relevant phenomenological 
role of the instanton contribution; as a further cross-check, we also show that S — > PP decays are badly described 
under a pure qq picture of the light scalar mesons, with or without instantons. The results are summarized in the 
Conclusions. 



2 Instanton effects in scalar meson dynamics 

QCD instantons produce an effective interaction which reduces the XJ(Nf)L x XJ(Nf)R global symmetry of the 
quark model in the chiral limit to SU (N/)l X S\J(Nf)n times baryon number. The effect can be described by the 
following effective Lagrangian [12] (see also [15]): 

Ci ex Det(Q LR ) , (Q LB r = qWr , (2) 

where i and j denote flavor indices (summation over color indices is understood) . 

With three light quark flavors, C\ is proportional to the product of three quark and three antiquark fields, 
antisymmctrised in flavor and color, and it includes a term of the type 

T r (jMj' 2 «l) , (3) 

where 

j\f = m<Mi, 4? ] = wi, (4) 

1 The fact that the J PC = ++ components of the qq states must have L = 1 can most easily be understood as follows. The operator 
q(0)q(0) can be written as q~LlR + QrQl- Now <?l an d qR both have helicity R but opposite momentum, hence opposite spin. Taking 
the momenta in the z-direction, we find that the spin structure of the fields here is (1, 2) + (2, 1), so 5* ot = 0. This is to be compared 
with the pion case, where 75 flips the relative sign: q~LQR — QrQl = (1, 2) — (2, 1). The latter is clearly the mode with total spin 5 = 0, 
since it can be written as eijUiVj, but in the scalar case, it is the S z = component of the state with total spin 5 = 1. To obtain 
total angular momentum J = we therefore need orbital angular momentum L = 1. The gamma algebra then shows that such L = 1 
contribution of the wave functions to <j(0)g(0) originates from the relativistic terms that mix the large components of the Dirac spinors 
with the small ones. 
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Fi gure 1'. The two main effects of the instanton Lagrangian in the scalar sector: (a) the tetraquark-^g mixing; (b) the Zweig-rule 
violating S — » PP amplitude. 



and 



[qq]ia = Ujkta^qflaq 1 ■ (5) 

is the spin-0 diquark operator, latin indices indicate flavor, greek indices stand for color and q c is the charge- 
conjugate of the quark field. As shown in Fig. [TJ this term can lead to a mixing between tetraquark and qq scalar 
states (Fig.[lJ,), and an effective S'(tetraquark)^ PP coupling (Fig. [TJd) which allows the ideally-mixed fo state in 
(fi]) to decay into two pions. 



2.1 The tetraquark— qq mixing 

The quantum numbers of j]^ 9 ' and jff^ match those of the tetraquark and qq scalar states $y and S: 



S' 



a |_ Jo 

a~ 

K~ 

a' | </"" 
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(S=[qq][qq]) , 



(S* = qq) , 



(G) 



(7) 



where the the neutral isoscalar states are not necessarily mass eigenstatcs. The instanton Lagrangian thus generate 
a mixing term 

£mix = 7 Tr(S • 5') . (8) 

A mixing of this form was introduced in [8] and is essential for a consistent identification of the scalar mesons 
around 1.5 GeV in terms of qq states, with a possible addition of a scalar glueball [17]. Indeed, the well identified 
1=1 and 1=1/2 states around 1.5 GeV, ao(1450) and K* (1430), also show a reversed mass ordering, although 
smaller than in the case of the light states. The anomaly can be explained as a contamination of the inverse 
hierarchy of the light tetraquark states via £ m i X - The coefficient 7 was determined phenomenologically from a fit 
to the mass spectrum [8] to be: 



|7| ~ 0.6 GeV^ 



(9) 



With this value and the observed masses, the bare masses of the lightest qq scalars turn out to be slightly above 
1 GeV. The result goes well with the estimate obtained in [13] from a consistent description of pseudoscalar states 
(including the 77') and scalar qq states within a linear sigma model. The bare qq masses agree also with the natural 
ordering of P-wave states, that predicts ++ masses to be smaller than 1 ++ and 2 ++ masses [8,18]. 

It was observed in [18] that the value of the effective coupling in ((9]) is much larger than what expected by 
usual QCD interactions for such Zweig-rule violating effect. To obtain this mixing by usual QCD interactions, it 
is necessary to annihilate completely quarks and antiquarks in the initial state and to produce from vacuum those 
of the final state. A strongly suppressed transition [18], which instead is provided almost for free by the instanton 
Lagrangian. 
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Figure 2: Leading quark-flavor diagrams for the decays into two pscudoscalar mesons of tetraquark (a) and qq (b) scalar mesons. 



A proviso concerns the hadronic matrix element of . In the fully non-relativistic approximation, one would 
have: 

(0|?(0)?(0)|5')=*(0), (10) 

where \I/(0) is the non-relativistic wave-function in the origin, which vanishes for P-waves. However, for relativistic 
quark fields we get a non- vanishing result, proportional to v — p/E. For QCD, Coulomb like, bound states, v ^ as, 
and the P-wave nature of S' results only in a mild suppression. 

The instanton induced mixing could in principle be determined by the following matrix element 

(S\d\S') oc (S\Tr{J [4q] J [2q] )\S') « Tr ((S\ J [4q] \0)(0\ J [2q] \S')^j . (11) 

However, at present we do not have a reliable independent information on the matrix elements of and J^ 2qi 
between scalar states and vacuum. 



2.2 S -> PP decays 

The leading mechanisms describing the decays of tetraquark and qq scalar states into two pseudoscalar mesons are 
illustrated in Fig.[2j In the tetraquark case, the diagram in Fig. [2^, (denoted as quark rearrangement in [9]) explains 
the affinity of fo and ao to KK channels. A non vanishing /o — > tttt amplitude appears if the possible mixing of 
/q ' and crt ! in JT]) is taken into account. However, as we illustrate in a quantitative way in the following section, 
the mixing alone does not lead to a good fit of all S — > PP decays [8,9]. As already observed in [9], this problem 
indicates the presence of additional contributions to the S — > PP amplitudes, generated by a different dynamical 
mechanism. This mechanism can be traced back to the instanton amplitude in Fig. [TJd. 

In constructing effective Lagrangians for the S(S') — ► PP decays an important role is played by chiral symmetry. 
In the chiral limit the octet components of the light pscudoscalar mesons 



$ = 



V2 



V2 



TT 

K 



7T U I <lq 

V2 ^ v/2 

K° 



K+ 
K° 

Vs 



(12) 



can be identified with the Goldstone bosons of the spontaneous breaking of SU(3)l x SU(3)a into SU(3)l+^. 
Following the general formalism of Ref. [19] (see also [20]), the effective Lagrangian of lowest dimension allowed 
by chiral symmetry contributing to S — ► PP decays, with P restricted to the octet components, turns out to be 
composed by only two independent operators: 



and 



2 (S) = — Tr(S)Tr(u^) 



(13) 



where 



1^ = ivJdjJJu' = — 



2/F)d fl $ + 0{<S> 2 



U 



_ „2 _ ^i^/F 



(14) 



Being dictated only by chiral symmetry, an identical structure holds for the 5" 
states. 



PP transitions of the qq scalar 
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The relative strength of the effective couplings of these operators can be determined by the correspondence of 
their flavor structures with a given quark-flavor diagram. In the tetraquark case, the leading amplitude in Fig. [2^ 
contributes to both 0% and O2, generating the combination 0\ — \Oi [9]. Taking into account also contributions 
with the singlet pscudoscalar field, the effective operator generated is: 



O f (S) = 1 (S)--0 2 (S) 



O s (S) , 



F 2 

Os(S) = — 



-Tr(S Ufi )Tr(u») + -Tr(S , )Tr(u p )Tr(u /i ) 



(15) 



The instanton induced coupling in Fig.[T}3 has a completely different flavor structure. From the chiral realization 
of the currents in (J3]) , 



(PP\jW\0) = (u^)o^) + --- 
(PP\Tr(J^J^)\S) = Tr(Su^) 0( ^ 



(16) 



where dots denote higher-order terms in <E> and in the chiral expansion. It follows that instanton effects are encoded 
only by Cq («!?). Taking into account both the leading quark- rearrangement diagram and the instanton contribution, 
decays of scalar tetraquark states into pscudoscalar mesons should be described by the effective Lagrangian 



decays^) = CfOf(S) + ClO^S) 



(17) 



where we expect |c/| <C |c/| since that the instanton contribution is a subleading effect. 

Similarly to the case of the S-S' mixing, we are not able to evaluate the hadronic matrix element of J' 49 ' 
in (|16[) from first principles, therefore we cannot predict the value of the effective instanton coupling c/ in Q17p . 
However, an interesting crosscheck of the normalization of the instanton effective Lagrangians is obtained under the 
hypothesis that the leading contribution to the amplitude in Fig. [TJd is the S' pole term arising by the contraction 
of the diagrams in Fig. [1^, and Fig. [2b- 

The S' — > PP decays of the heavy qq states have been analyzed in Ref. [8] and found to be reasonably well 
described by the effective Lagrangian 

-Cdecay (SO = c' f O^S') , (18) 

corresponding to the chiral realization of the diagram in Fig. [2}d. The value of the effective coupling is found to be 
Cj fa 6.1 GcV -1 (cj = 2A' in the notation of Ref. [8]). Under the plausible hypothesis of S' pole dominance for 
the instanton-induced S — > PP amplitude we thus expect 



^(5'— pole) 



bl 



Ml 



l7l df 
M% 



0.016 



(19) 



where we have used Mg/ ~ 1.5 GeV. 

The above estimate should be taken cautiously. Non-pole terms are not expected to be totally negligible, and 
the poor experimental information about the heavy scalar meson decays implies a sizable uncertainty in the value 
of Cj. Nonetheless it is very encouraging that the value in (fT9|) is consistent with the phenomenological value 
of cj that we shall obtain in the next section from the phenomenology of light scalar decays with the effective 
Lagrangian l[T7|). 



3 Phenomenological analysis of S — ► PP decays 

In this section we analyze the decays of the light scalars into two pseudoscalar mesons using the effective Lagrangian 
in p7p . As discussed before, such Lagrangian with |c/| |c/| corresponds to the tetraquark hypothesis for the 
light scalar mesons, taking into account instanton effects. 
The g SPl p 2 couplings, defined by 

A(S - Pi( P i)P 2 (p 2 )) = .g SPlP2 = g SPlP2 \{Ml - M 2 Pi - Ml 2 ) , (20) 



5 





a 


K 


fo 


a 


M(MeV) 


441 


800 


965 


999 



Table 1: Numerical values used for the light scalar masses (sec Rcf. [1—4,22]) 



are reported in Table [2l The decay rates are then given by: 

r(5^PxP 2 ) = \A(S 



P1P2) 



8ttM| 



where p* is the decay momentum. 

For comparison, we also attempt a fit of light scalar meson decays with the effective Lagrangian 

£ cS (S') = c' f Ox(S') + c'jOfiS') (21) 

where, contrary to what advocated so far, we identify the qq scalar field in {7} with the nonet of light scalar mesons. 
In this scheme the leading operator is Oi(S') while instanton effects contribute to Of(S') (the corresponding g SPl p 2 
arc reported in the third line of Table [2|). The comparison shows that, beside the problems with the mass spectrum, 
the qq hypothesis for the light scalar mesons has also serious difficulties in fitting S — > PP data. 



3.1 Mass mixing of neutral states 

The physical /o and a states are in general a superposition of the ideally mixed states and /g ' defined in 
Eqs. ©-((T]). Introducing a generic mixing between the mass eigenstates and the ideally mixed states 



a 




cosw 


sin 10 


A 




— sin oj 


cos LJ 



Jo 



(22) 



the value of the mixing angle u is determined by the experimental values of the scalar masses (see [9,21]). Fixing 
the scalar masses to the input values in Table [T] but for the poorly known k mass, and letting the latter vary in 
the interval [750, 800] MeV, we find that the deviations from the ideal mixing case are quite small. We find in 
particular \lu\ < 5°, which we use as range to estimate the impact of a non- vanishing to. 

Another important ingredient to compute the physical amplitudes is the 77-77' mixing. In the octet-singlet basis 
we define 



= U{4> PS ) 



U(4> P3 ) 



COS( 

sin < 



PS 
'PS 



— sin 

cosd 



rps 
PS 



It is useful to consider also the mixing in the quark basis: 



n 



U{-6) 



(23) 



(24) 



where <p ps + 9 + tan 1 (v / 2) = 0. From the analysis of the pscudoscalar meson masses, 77 decays of r\ and 77' and 
J/ip -> 7 77/77' [26], one obtains <j) PS ~ -22° (0 ~ -33°). 

The effect of 77 — 77' mixing in the ao — > tjtt processes plays an important role in distinguishing the two hypotheses 
for the scalar mesons. In the four-quark case, the quark exchange amplitude (Fig. [2jt) produces a pure 77.,, while 
the instanton interaction (Fig. [TJs) produces a pure rj q . Expressing the coupling to physical particles 77 and 77' in 
terms of the octct-singlct mixing angle leads to: 



AH 



= V2, 



ci cos0 — Cf smf 



v2(c/ + ci) cos (j) ps + (cf - 2c/)sin</> j: 

7s 



In the qq hypothesis, where the role of two operators 0\ and Of is exchanged, one finds 

y/2(df + c'j) cos(/) ps - (2df - c'j) sin0 pg 



9 [ f + = 

J dg 7T+J) 



V3 



(25) 



(26) 
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c' f + C'j 


V2c> 


-di 


c' f + c'j 



Tabic 2: The g SPlP ^ couplings for tetraquark and qq scalar mesons, from the effective Lagrangians 1 11 71 1 and (1211 . respectively. 



Processes 


KLOE [23,25] 


BES, Crystal Barrel, WA102 


fo -> X+X- 


i.43i2;22 (i.3 ±o.i) 

3.761^ (0.4lg;S) 


2.32 ±0.25 [27] 
4.12 ±0.55 [27] 


a — > 7r°?7 
a -> K+if- 


2.8 ±0.1 (2.2 ±0.1) 
2.16 ±0.04 (1.6 ±0.1) 


2.3 ±0.1 [28] 2.1 ±0.23 [29] 
1.6 ±0.3 [28] 



Table 3: Experimental data for the S — > PP amplitudes in GeV. The number in brackets in the KLOE column refers to the 
parameterization of Ref. [30] without the a pole (the absence of the a contribution makes the fo — > 7r+7r — , K + K~ results between 
brackets less reliable for the present analysis; similarly, we do not quote the fo — > 7r+7r~ , K + K~ results extracted from <r(e + e - — > 
7r + 7r _ 7) [24] which suffer of a larger background). The numbers in the last column refer to other experiments where it has been possible 
to unambiguously extract the information on the partial amplitudes. 



3.2 Numerical analysis 

The results of the fit to S — > PP amplitudes are reported in Table 2J We use the masses given in Table Q] and 
<f) ps = —22° [26]. We analyze both the two- and four-quark hypotheses, with or without the instanton contribution, 
using a — > tttt and fo — > tttt as input channels. For simplicity, in Table [4] we compare the theoretical predictions 
with the KLOE data only. A comparison with other experimental results can simply be obtained by means of 
Table [3l where all the available experimental information is collected. The data for a and K decays are the results 
of the theoretical analyses in [3,4]. 

The values of the couplings derived in the four-quark hypothesis and w = are: 

c f = 0.041 MeV" 1 ; a = -0.0022 MeV" 1 (cj = 0) (27) 

The negative sign of cj is chosen to minimize the ao — > t/tt rate. As shown by eq. (|25j) . the negative interference 
between c/ and cj increases for (p ps more negative. 

As an alternative strategy, we have performed a global fit, assigning conventionally a 10% error to the a rate, 
15% to k, 30% to all others and <fi PS = —22° and searched for a best fit solution. The results of this fit are reported 
in the fourth column. The best-fit couplings are: 

c f = 0.020 ±0.002 McV" 1 , c/ = -0.0025 ± 0.0012 MeV" 1 (-5° < u < 5°) (28) 

Central values are for u> = 0, errors are given by letting to to vary in ±5° range (see sect. 13.11) . 
Some comments concerning the fit under the four-quark hypothesis arc in order: 

• There is a good overall consistency, the fit is stable and, as expected, |c/| -C \cf\. The fitted value of |cj| is 
also perfectly consistent with the pole estimate presented in Eq. (fT9|) . 

• The positive feature of the instanton contribution is that it provides a non-vanishing foTnr coupling inde- 
pendent from mixing and, at the same time, it improves the agreement with data on the 'clean' ao — ► rjir 
channel. 

• The relation between k — > Kir and ao — > KK is fixed by SU(3) and does not depend on the value of the 
couplings. As a result it is impossible to fit simultaneously the central values of these two amplitudes without 
introducing symmetry breaking terms. 
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Processes 


Ah([<7<7] [<?<?]) 




Ah te?) 






with inst. 


no inst. 


best fit 


with inst. 


no inst. 




a — > 7T + 7T~ 


input 


input 


1.6 


input 


input 


3.22 ±0.04 


K+ -> K\+ 


7.3 


7.7 


3.3 


6.0 


5.5 


5.2 ±0.1 


fo -> TT+TT- 


input 


[0-1.6] 


1.6 


input 


[0-1.6] 


1.4 ±0.6 


fo -> X+X- 


6.7 


6.4 


3.5 


6.4 


6.4 


3.8 ± 1.1 


a — ► 7T U 77 


6.7 


7.6 


2.7 


12.4 


11.8 


2.8 ±0.1 


a -► if+if" 


4.9 


5.2 


2.2 


4.1 


3.7 


2.16 ±0.04 



Table 4: Numerical results, amplitudes in GeV. Second and third columns: results obtained with a decay Lagrangian including or not 
including instanton effects, respectively. Fourth column: best fit with instanton effects included (see text). Fifth and sixth columns: 
predictions for a qq picture of the light scalars with and without instanton contributions. All results are obtained with a r]—r)' mixing 
angle <p ps = —22°. Second and fifth columns are computed with a scalar mixing angle w = 0. The fo — > tttt couplings in the third 
and sixth columns are computed with u in the range ±5° (see text). Data for cr and k decays are from [3,4], the reported amplitudes 
correspond to: r t ot(<r) = 544 ± 12, r t ot(«) = 557 ± 24. 



• Being extracted from cr(e + e — > tttt^) data, the experimental values for the fo, ao — > KK couplings reported 
in Table [3] are subject to a sizable theoretical uncertainty (not shown in the table). 

• A mixing angle u in the ±5° range (see sect. I3.1[) affects, most significantly, the no-instanton case, leading 
to a non-vanishing fo — >■ tttt amplitude marginally consistent with data but it does not solve the ao — ► rjTT 
problem (Table [4j third column). 

The fit under the qq hypothesis, reported in the fifth and sixth columns of Table |4] is much worse. Due to the 
exchange of fo and a going from S to S', it remains true that there is no fo — ► tttt in the absence of instant ons 
and, similarly to the tetraquark case, we obtain \c\\ <C \dA. However, the situation is drastically different for the 
ao —> TT°ri channel. In the qq case there is no sign conspiracy that produces the cancellation found before and the 
predicted amplitude is far from the observed one for any value of <j> PS , with or without instanton effects. This bad 
fit provides a further evidence, beside the inverse mass spectrum, against the qq hypothesis for the lightest scalar 
mesons. 

4 Conclusions 

The addition of the instanton-induced effective six-fermion Lagrangian lead us to a simple and satisfactory descrip- 
tion of both the light scalar mesons below 1 GeV and the heavier scalar states around 1.5 GeV. The light mesons 
are predominantly tetraquark states (£), while the heavier ones are predominantly qq states (£')■ Instantons in- 
duce a mixing between the two sets of states, which explains the puzzling mass spectrum of the heavier mesons. 
Moreover, integrating out the heavy states, instanton effects manifest themselves in the dynamics of the lightest 
states, generating the Zweing-rule violating amplitude which is necessary for a consistent description of the strong 
decays of the light scalar mesons. 

The phenomcnological determinations of the S-S' mixing and of the Zwcig-rulc violating S — > PP amplitude 
suggest that the dynamics of the scalar mesons could be described, to a good extent, by a simple effective Lagrangian 
of the type 

£eff,aii = Tr (SM 2 S S) + Tr (S'M%,S') + 7 Tr (SS') + c s O s (S) + c'p^S') (29) 

where M. 2 S{I ) = + are appropriate mass matrices reflecting the inverse and normal mass ordering of 

tetraquark and qq states, and 7 is the coupling encoding the instanton contribution. 

Note Added After this work was submitted, an independent analysis of the instanton-induced mixing between 
two and four quark scalar mesons has been presented in [31]. 
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